Signature of Gravity Waves in Polarization of the Microwave Background 
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Using spin-weighted decomposition of polarization in the Cosmic Microwave Background (CMB) 
we show that a particular combination of Stokes Q and U parameters vanishes for primordial fluc- 
tuations generated by scalar modes, but does not for those generated by primordial gravity waves. 
Because of this gravity wave detection is not limited by cosmic variance as in the case of tempera- 
ture fluctuations. We present the exact expressions for various polarization power spectra, which are 
valid on any scale. Numerical evaluation in inflation-based models shows that the expected signal 
is of the order of 0.5 fJi,K, which could be directly tested in future CMB experiments. 
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It is now well established that temperature anisotropies 
in CMB offer one of the best probes of early universe, 
which could potentially lead to a precise determination 
of a large number of cosmological parameters . The 
main advantage of CMB versus more local probes of 
large-scale structure is that the fluctuations were created 
at an epoch when the universe was still in a linear regime. 
While this fact has long been emphasized for temperature 
anisotropies, the same holds also for polarization in CMB 
and as such it offers the same advantages as the temper- 
ature anisotropies in the determination of cosmological 
parameters. The main limitation of polarization is that 
it is predicted to be small: theoretical calculations show 
that CMB will be polarized at 5-10% level on small angu- 
lar scales and much less than that on large angular scales 
[|,|. However, future CMB missions (MAP, Planck) will 
be so sensitive that even such low signals will be measur- 
able. Even if polarization by itself cannot compete with 
the temperature anisotropies, a combination of the two 
could result in a much more accurate determination of 
certain cosmological parameters, in particular those that 
are limited by a finite number of multipoles in the sky 
(i.e. cosmic variance). 

Primordial gravity waves produce fluctuations in the 
tensor component of the metric, which could result in a 
significant contribution to the CMB anisotropies on large 
angular scales. Unfortunately, the presence of scalar 
modes prevents one from clearly separating one contri- 
bution from another. If there are only a finite number of 
multipoles where tensor contribution is significant then 
there is a limit in amplitude beyond which tensors can- 
not be distingushed from random fluctuations. In a noise 
free experiment the tensor to scalar ratio T / S needs to be 
larger than 0.15 to be measurable in temperature maps 
1^. Independent determination of the tensor spectral 
slope riT is even less accurate and a rejection of the con- 
sistency relation in inflationary models T/S ~ —Iut is 
only possible if Itt-tI ^ {T/S)/7 ||,^. Polarization pro- 
duced by tensor modes has also been studied pi , but only 



in the small scale limit. In previous work correlations be- 
tween Stokes parameters Q and U have been used. These 
two variables are not the most suitable for the analysis 
as they depend on the orientation of coordinate system. 
It was recently shown that in Fourier space Q and 
U can be decomposed in two components, which do not 
depend on orientation. Moreover, scalar modes only con- 
tribute to one of the two, leaving the other as a probe of 
gravity waves. These arguments have been made in the 
small angle approximation. In this Letter we remove this 
limitation by presenting a full spherical analysis of polar- 
ization using Newman-Penrose spin-s spherical harmonic 
decomposition. An alternative decomposition in terms 
of tensor harmonics has been presented recently by 
We show that there is a particular combination of Stokes 
parameters that vanishes in the case of scalar modes, 
which can thus be used as a probe of gravity waves. We 
present the expression for the power spectrum of various 
polarization components using the integral solution |^ 
and evaluate it numerically for a variety of cosmological 
models. We also discuss the sensitivity needed to detect 
this signal and compare it to the expected sensitivities of 
future CMB satellites. 

Linear polarization is a symmetric and traceless 2x2 
tensor |l^ that requires 2 parameters to fully describe it: 
Q and U Stokes parameters. These parameters depend 
on the orientation of the coordinate system on the sky. 
It is convenient to use Q + iU and Q — iU as the two 
independent combinations, which transform under right- 
handed rotation by an angle 4> a,s {Q + iU)' = e~^*"^((3 -I- 
iU) and (Q - iU)' = e2*"^(Q - iU). These two quantities 
therefore have spin- weights 2 and —2 respectively and can 
be decomposed into spin ±2 spherical harmonics ±2i^m 
(for a discussion of spin- weighted harmonics see [pd|) 

(Q + zC/)(n) = ^ 

hn 

{Q - iU){fl) = ^ a-2,lm-2Ylm{.'h). (1) 
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Spin s spherical harmonics form a complete orthonormal 
system for each value of s. An important property of 
spin-weighted basis is that there exists spin raising and 
lowering operators & and 6 (see for their explicit 
form). By acting twice with a spin lowering and rais- 
ing operator on {Q + ill) and {Q — ill) respectively one 
obtains quantities of spin 0, which are rotationally invari- 
ant. These quantities can be treated like the temperature 
and no ambiguities connected with the orientation of co- 
ordinate system on the sky will arise. Conversely, by 
acting with spin lowering and raising operators on usual 
harmonics spin s harmonics can be written explicitly in 
terms of derivatives of the usual spherical harmonics . 
Their action on ±2Yim leads to 



overall sign to agree with the small scale expressions in 
M . Note that our E and B are proportional to G and C 

rTHi (S) (S) 

We find that — and only ai^ is non-zero. 
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The polarization power spectrum is defined as the rota- 
tionally invariant quantity C; = 2lTT "^^"^ 
E its ensemble average can be obtained by acting twice 
with spin raising (or lowering) operator on equation ^ 
leading to (see llj] for details) 
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6'{Q-^U){n) = J2( K±|l:y^'a_2,,„,li„,(n). (2) 



where ji{x) is the spherical Bessel function of order I 
and Pfjiik) is the primordial power spectrum of scalar 
metric perturbations, usually assumed to be a power law 
P^{k) (X fc"--*. 

In the case of tensors the form for Q and U in the 
frame where A; || 2; is 



With these definitions the expressions for the expansion 
coefficients of the two polarization variables become 



Q{n,k) ^ -{1 + fi'^)e'^'^ J dTe''^''g'i{k,T) 
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dn Y;^in)6'iQ + zU)in) 
dnY;^{n)6'{Q-iU){n). (3) 



Uin,k) = -2ifie'"l' J dTe"''c/*(fc, t) + c.c, (6) 

where the source is a complex sum over the two indepen- 
dent tensor polarization states, 4" — (5"+ — i5'^)/2, and 
can be expressed in terms of temperature and polariza- 



To obtain the expression for the polarization power 
spectrum we will use the integral solution of the Boltz- 
mann equation In the case of scalar perturbations 
for any given Fourier mode k only Q^"^' is generated in 
the frame where k \\ z |l2[. 
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6 '^{Q — ill) are not equal, so both and a 

be nonzero. Using a similar procedure as above we obtain 
their power spectra fill] 



Q(^)(n, k) = J(l - ^,') I dre"''5(r)n(fc, r) (4) C = {Att)' J k'dkP^{k) 



where x — k{TQ — t) and r is the conformal time with 
To its present value. Directions in the sky are denoted 
with polar coordinates {9, cf)) and n = cos{d). We in- 
troduced the visibility function ^(t) = ke~'^\ where k 
is the differential optical depth for Thomson scatter- 
ing, k = aUeXeCTT^ a(T) is the expansion factor nor- 
malized to unity today, rig is the electron density, Xe 
is the ionization fraction and ax is the Thomson cross 
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section. The source term 11 = + A 
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ApQ was 



expressed in terms of temperature quadrupole , po- 

larization monopole ApQ and its quadrupole Ap2 . Be- 
cause U'-^^ ~ and Q^'^-' is only a function of /i in the 
k II z frame it follows S^{Q + iU) = 6 ^{Q - iU) and 
so a2^i^ — ;m- convenient to introduce two or- 
thogonal combinations aE,im = ~{a2dm + a-2,im)/2 and 
a-B.im = (a2,Zm— a-2,im)/2. Here i? and _B refer to electric 
and magnetic type parities ||l^] and we have chosen the 



where Ph{k) oc fc"^^'^ is the primordial power spectrum 
of gravity waves. In the small scale limit these expres- 
sions agree with those derived previously . 

Using the above expressions we may numerically eval- 
uate the power spectra in various theoretical models. We 
use T/S as the parameter determining the amplitude 
of tensor polarization. Fig. |l| shows the predictions for 
scalar and tensor contribution in standard CDM model 
with no reionization (a) and in reionized universe with 
optical depth of k = 0.2 (b). The latter value is typical in 
standard cosmological models 111] . We assumed T/S — 1 
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and riT = {ng — 1) = —0.15. In the no-reionization case 
both tensor spectra peak around I ~ 100 and give com- 
parable contributions, although the B channel is some- 
what smaller. Comparing the scalar and tensor E chan- 
nels one can see that scalar polarization dominates for 
T/Si 1. Even though tensor contribution is larger than 
scalar at low I, the overall power there is too small to be 
measurable. Tensor reconstruction in the E channel suf- 
fers from similar drawbacks as in the case of temperature 
anisotropics: because of large scalar contribution cosmic 
variance prevents one to isolate very small tensor contri- 
butions 1^ . The situation improves if the epoch of reion- 
ization occured sufficiently early that a moderate optical 
depth to Thomson scattering is accumulated (Fig |l|b). 
In this case there is an additional peak at low I ||l^ and 
the relative contribution of tensor to scalar polarization 
in E channel around Z = 10 is higher than around / — 100. 
Still, if T/S <C 1 cosmic variance again limits one to ex- 
tract unambiguously the tensor contribution. It is in this 
limit that the importance of B channel becomes crucial. 
This channel is not contaminated by scalar contribution 
and is only limited by noise, so in principle with sufficient 
noise sensitivity one can detect even very small tensor to 
scalar ratios. Moreover, a detection of signal in this chan- 
nel would be a model independent detection of non-scalar 
perturbations. In the following we will discuss sensitivity 
to gravity waves using both only B channel information 
and all available information. 

We can obtain an estimate of how well can tensor pa- 
rameters be reconstructed by using only the B chan- 
nel and assuming that the rest of cosmological param- 
eters will be accurately determined from the tempera- 
ture and E polarization measurements. While this test 
might not be the most powerful it is the least model 
dependent: any detection in B channel would imply a 
presence of non-scalar fluctuations and therefore give a 
significant constraint on cosmological models. Because 
the B channel does not cross-correlate with either T or 
E only its auto-correlation needs to be consid- 

ered. A useful method to estimate parameter sensitivity 
for a given experiment is to use the Fisher information 
matrix §,1,1,0 
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where /sky is the sky coverage. Receiver noise can be 
parametrized by A'kg'^ /N , where a is the noise per pixel 
and N is the number of pixels. Typical values are 
(0.15/xK)2 for MAP and {Qmb^Kf for the most sensi- 
tive Planck bolometer channel in one year of observation. 
In our case the parameters Si can be T / S and ut, so that 
the matrix is only 2x2. The error on each parameter is 



given by {cx~^Y^'^ if the other parameter is assumed to be 
unknown and {aii)~^/'^ if the other parameter is assumed 
to be known. Using this expression we may calculate the 
experiment sensitivity to these parameters. Current in- 
flationary models and limits from large scale structure 
and COBE predict T/S to be less than unity. Figure |l| 
shows that the expected amplitude in this case is below 
0.5/iK. We find that MAP is not sufRciently sensitive 
in B channel to detect these low levels. On the other 
hand, Planck will be much more sensitive and can detect 
T/S > 0.3 if tensor index ut is assumed to be known (for 
example through the consistency relation). For the un- 
derlying model with T/S = 1 one can determine it with 
an error A{T/ S) ^0.1. If tensor index is not known then 
a combination of the two parameters, which corresponds 
to the total power under the B curve in figure can still 
be determined with the same accuracy. 
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FIG. 1. Multipole moments for the three polarization spec- 
tra for no-reionization case (a) and reionized case with optical 
depth of 0.2 (b). The underlying model is "standard CDM" 
with T/S = 1. 

Separate determination of the tensor amplitude and 
slope from the B channel is only possible in reionized 
models. In the no-reionization model the contribution to 
B is very narrow in / space and the leverage on ut inde- 
pendent of T/ S is small, so that the correlation coefficient 
ai2/(cnia22)^^^ is almost always close to unity. A mod- 
est amount of reionization improves the separation; in the 
reionized models the power spectrum for B is bimodal 
(figure |l|) and the overall signal is higher, which gives a 



3 



better leverage on tit independent of T/S. For k = 0.2 
the Planck errors are A(r/5) ^ 0.15 and Any ~ 0.1 for 
the underlying model with T/S — 1. These results de- 
pend on the overall amplitude relative to the noise level. 
As long as both peaks can be separated from the noise 
one can determine the tensor slope, which allows to test 
the inflationary consistency relation. 

Combining all the information by adding temperature, 
E polarization and their cross-correlation further im- 
proves these estimates. In this case other parameters 
that affect scalar modes such as baryon density, Hubble 
constant or cosmological constant enter as well and the 
results become more model dependent [Q. Fisher infor- 
mation matrix has to be generalized to include all the 
parameters that can be degenerate with the tensor pa- 
rameters. The results depend on the class of models and 
number of parameters one restricts to in the analysis, as 
opposed to the results based on B channel above, which 
depend only on the two main parameters that charac- 
terize the gravity wave production. As a typical ex- 
ample, for T/S = 0.1 and k = 0.1 one can determine 
A{T/S) = 0.05 and Aut = 0.2 with Planck §. These 
errors improve further if a model with higher T/S or n 
is assumed. For the same underlying model without us- 
ing polarization the expected errors are AT/S ^ 0.26 
and Aut ~ 1, significantly worse than with polarization. 
Even for MAP the limits on T/S improve by a factor of 
2 when polarization information is included. 

To summarize the above discussion, future CMB mis- 
sions are likely to reach the sensitivities needed to mea- 
sure (or reject) a significant production of primordial 
gravity waves in the early universe through polarization 
measurements, which will vastly improve the limits from 
temperature measurements only and allow a test of con- 
sistency relation. The more challenging question is the 
foreground subtraction at the required level. At low fre- 
quencies radio point sources and synchrotron emission 
from our galaxy dominate the foregrounds and both are 
polarized at a 10% level. Their contribution decreases at 
higher frequencies and with several frequency measure- 
ments one can subtract these foregrounds at frequencies 
around 100 GHz at the required microkelvin level. At 
even higher frequencies dust is the dominant foreground, 
but is measured to be only a few percent polarized . 
We hope that the signature of gravity waves discussed 
here would provide further motivation to pursue the fea- 
sibility studies of polarization measurements. 

While we only discussed scalar and tensor modes, vec- 
tor modes, if present before recombination, will also con- 
tribute to both polarization channels and so could con- 
taminate the signature of gravity waves. At present there 
are no viable cosmological models that would produce a 
significant contribution of vector modes without a com- 
parable amount of tensor modes. In inflationary models 
vector modes, even if produced during inflation, decay 
away and are not significant during recombination. In 



topological defect models nonlinear sources continuously 
create both vector and tensor modes and so some of the 
signal in B channel could be caused by vector modes. 
Even in these models however some fraction of signal in 
B will still be generated by tensor modes and in any 
case, absence of signal in B channel would rule out such 
models. Polarization thus offers a unique way to probe 
cosmological models that is within reach of the next gen- 
eration of CMB experiments. 
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